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The 16 0( 3 He,d) 17 F reaction has been used to determine asymptotic normalization coefficients for 
transitions to the ground and first excited states of 17 F. The coefficients provide the normalization 
for the tails of the overlap functions for 17 F — > ie O + p and allow us to calculate the S-factors 
for 16 0(p,^) 17 F at astrophysical energies. The calculated S-factors are compared to measurements 
and found to be in very good agreement. This provides the first test of this indirect method to 
determine astrophysical direct capture rates using transfer reactions. In addition, our results yield 
S(0) for capture to the ground and first excited states in 17 F, without the uncertainty associated 
with extrapolation from higher energies. 
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I. INTRODUCTION 

Nuclear capture reactions, such as (p, 7) and (01,7), 
play a major role in defining our universe. A primary goal 
in nuclear astrophysics is to determine rates for capture 
reactions that are important in the evolution of stellar 
systems. However, the reactions of interest often involve 
radioactive targets which makes measurements quite dif- 
ficult or even impossible using conventional methods. 
Hence techniques have been developed to determine rates 
by indirect methods. For example, precise information on 
excitation energies and particle decay widths can be used 
to make accurate predictions of rates which proceed by 
resonance capture. The only reliable method to deter- 
mine a reaction rate that is dominated by direct capture 
has been to measure it at laboratory energies with a low 
energy particle beam and then extrapolate the result to 
energies of astrophysical interest. 

Attempts have been made to use both Coulomb dis- 
sociation jl| and the determination of asymptotic nor- 
malization coefficients (ANC) from conventional nuclear 
transfer reactions | [|] to determine S-factors for direct 
capture reactions, but neither technique has been tested 
to verify its reliability. Such tests are crucial, as stressed 
in the most recent evaluation of solar fusion cross section 
rates We report here the first test of one of these 
two techniques to determine astrophysical S-factors; we 
demonstrate that the ANC inferred from a measurement 
of a proton transfer reaction can directly determine a 
(p, 7) direct capture rate at astrophysical energies. 

Direct capture reactions of astrophysical interest usu- 
ally involve systems where the binding energy of the cap- 
tured proton is low. Hence at stellar energies, the cap- 
ture proceeds through the tail of the nuclear overlap wave 
function. The shape of this tail is completely determined 



by the Coulomb interaction, so the rate of the capture 
reaction can be calculated accurately if one knows its 
amplitude. The asymptotic normalization coefficient C 
for the system B <-> A + p specifies the amplitude of 
the single-proton tail of the wave function for nucleus 
B when the core A and the proton are separated by a 
distance large compared to the nuclear radius. Thus, 
this normalization coefficient determines the correspond- 
ing direct capture rate. 

The advantage of the ANC approach is that it pro- 
vides a method to determine direct capture S-factors ac- 
curately from the results of nuclear reactions such as 
peripheral nucleon transfer which can be studied with 
radioactive beams and have cross sections that are or- 
ders of magnitude larger than the direct capture reac- 
tions themselves. Furthermore, direct capture S-factors 
derived with this technique are most reliable at the lowest 
incident energies, precisely where capture cross sections 
are smallest and most difficult to measure directly. In 
fact, the ANC approach even permits one to determine 
S-factors at zero energy, which is not possible with direct 
measurements except by extrapolation. 

While there is little controversy that knowledge of the 
asymptotic normalization coefficient for a loosely bound 
nuclear system allows one to compute the correspond- 
ing direct capture rate, the nuclear astrophysics commu- 
nity has clearly indicated || that a test of the relation- 
ship between the transfer reaction cross section and the 
astrophysical S-factor is important to validate this ap- 
proach. The community's skepticism originates in the 
well-known model dependence found in distorted-wave 
Born approximation (DWBA) analyses of transfer reac- 
tion data in terms of spectroscopic factors, which is due 
to the uncertainty in the DWBA calculations associated 
with the choice of optical model potentials and single 
particle wave functions. By parametrizing the DWBA 
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cross section of a peripheral transfer reaction in terms of 
ANC's, rather than spectroscopic factors, we can reduce 
the uncertainty associated with the choice of single parti- 
cle wave functions so that it becomes small compared to 
that associated with the optical potential . By choos- 
ing appropriate reactions, beam energies and scattering 
angles, we can also minimize the uncertainty associated 
with the choice of optical model potentials. 

In this article, we describe a measurement of the 
w O( 3 He, d) 17 F reaction, from which we determine the 

ANC's for the | ground state and the \ first ex- 
cited state in 17 F. We then use our measured ANC's 
to calculate, with no additional normalization factors, 
the S-factors for the w O(p, 7) 17 -? 7, reaction at astrophys- 
ical energies. Such a determination of the S-factors for 
16 0(p, 7) 17 F from its ANC's measured in proton transfer 
reactions is an ideal test case for this indirect method || 
because the results can be compared to existing direct 
measurements of the capture cross sections |9|,|l0|]. Fur- 
thermore, the 16 0(p,j) 17 F reaction has substantial sim- 
ilarities to the 7 Be(p, j) 8 B reaction, which is the source 
of all high energy neutrinos produced in the sun. We 
will report determinations of the S-factor for 7 Be(p, 7) 8 i? 
using this technique in future publications. It will also 
be straightforward to utilize this procedure to determine 
S-factors at astrophysical energies for other cases that 
include significant direct capture components. 

II. 17 F ^ ie '0 + p ASYMPTOTIC 
NORMALIZATION COEFFICIENTS 

For a peripheral transfer reaction, ANC's are extracted 
from the measured angular distribution by comparison 
to a DWBA calculation. Consider the proton transfer 
reaction a + A — > c + B, where a = c + p and B = A+p. 
The experimental cross section is related to the DWBA 
calculation according to 

iBjBlaja 

where 

p _ <T l B jBlgja fi)S 

n-lBjBlaja ~ l2 U.2 ' V Z 7 

Apl B jB cpl a j a 

a is the calculated DWBA cross section and the b's are 
the asymptotic normalization constants for the single 
particle orbitals used in the DWBA. The sum in Eq. (|l|) is 
taken over the allowed angular momentum couplings, and 
the C's are the ANC's for B -> A+p and a -> c+p. The 
normalization of the DWBA cross section by the ANC's 
for the single particle orbitals makes the extraction of the 
ANC for B — > A + p insensitive to the parameters used 
in the single particle potential wells 0]q], in contrast to 
traditional spectroscopic factors. See jpfl for additional 
details. 



DWBA calculations of the 16 0( 3 He, d) 17 F reaction 
populating the 17 F first excited state indicate that the 
sensitivity of the extracted ANC to the choice of opti- 
cal model potentials is minimized near 0°. There ex- 
ists a previous study of the le O( 3 He,d) 17 F reaction at 
E3 He = 25 MeV |T| that reported cross sections at 9 
angles over the range 9 cm w 6 — 36°. The limited small- 
angle coverage makes an attempt to infer the 17 F first 
excited state ANC from that experiment very imprecise. 
We have now measured the le O( 3 He, d) 17 F reaction at 
Ea He ps 29.7 MeV primarily to determine the angular dis- 
tribution carefully at small angles, thus minimizing the 
systematic uncertainty in the extracted ANC. However, 
by obtaining data at a second beam energy, we can also 
do a combined analysis to reduce our sensitivity to the 
choice of optical potentials even further. 

Two separate measurements were performed, one op- 
timized to determine the absolute cross section with a 
minimum of uncertainty and the other to obtain a de- 
tailed angular distribution at small angles. The reaction 
was measured at laboratory angles between 6.5° and 25° 
using a momentum-analyzed 29.75 MeV 3 He beam from 
the U-120M isochronous cyclotron of the Nuclear Physics 
Institute (NPI) of the Czech Academy of Sciences inci- 
dent on a Mylar target. The target thickness was mea- 
sured to be 134 fig/ cm 2 by scanning with well-collimated 
alpha particles from 241 Am, 238 Pu and 244 Cm. Reac- 
tion products were observed by a pair of detector tele- 
scopes, consisting of 150 /im thick AE and 2000 fim thick 
E Si surface barrier detectors, with solid angles of 0.23 
msr. One of the telescopes was rotated about the target 
during the measurements while the other was fixed at 
6l = 18.2°. Elastic scattering and several reaction chan- 
nels were measured simultaneously in both telescopes to 
provide a continuous calibration of the beam energy, re- 
action angle and target thickness. The beam current was 
integrated by a Faraday cup biased to 1 kV. Absolute 
cross sections were determined to ±4.5%, using proce- 
dures developed at NPI to minimize overall normaliza- 
tion uncertainties |lj,[l3|]. 

Small angle data at laboratory angles between 1° and 
11° were obtained using a molecular ( 3 He — d) + beam 
from the Texas A&M University K500 superconducting 
cyclotron incident on a 540 /ig/cm 2 Mylar target. The 
angular spread of the beam on target was «0.1° after 
passing through the Texas A&M Beam Analysis System 
PU . Reaction products were detected at the focal plane 
of the Multipolc Dipolc Multipole magnetic spectrome- 
ter jyl using the modified Oxford detector j^6|. The 
detector consists of a 50 cm long gas ionization cham- 
ber to measure the specific energy loss of particles in 
the gas and their focal plane positions at four resistive 
wires, separated by 16 cm steps along the particles' tra- 
jectories, followed by an NE102A plastic scintillator to 
measure the residual energy. The 3 He energy in the 
molecular beam was determined from the crossover be- 
tween the 12 C( 3 He, t) 12 N and w O( 3 He, a) 15 reactions, 
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FIG. 1. Angular distributions for the ground and first ex- 
cited states of 17 F from the lfi O( 3 He, d) 17 F reaction. The 
dashed curves are DWBA fits using optical potential set I in 
Table I, and the solid curves use set II. 



observed simultaneously off the Mylar target. It was 
29.71 MeV, tuned to match the measurements carried 
out at the NPI. The beam angle was determined to ±0.1° 
from the crossover between the 1 H( 3 He 3 He) 1 !! and 
12 C( 3 iTe, 3 He) 12 C*(A.U MeV) reactions, also observed 
simultaneously off the Mylar target. The charge in the 
beam was collected in a Faraday cup and provided the 
normalization between different scattering angles. The 
spectrometer has an acceptance of A8l — 4°, which was 
divided into 8 separate 0.5° angle bins by ray tracing. It 
was moved in 2° steps from laboratory angles of 3° to 
9°. With this procedure, the internal consistency of the 
normalization between angles was verified. Additional 
details regarding the experimental procedures may be 
found in . 

The absolute normalization of the Texas A&M cross 
section measurements was determined by matching the 
ground and first excited state yields to those determined 
at NPI in the angular region where the two data sets 
overlap. The matching procedure introduced an addi- 
tional ±1.1% uncertainty in the absolute normalization 
of the small angle cross section measurements. The com- 
bined angular distributions for the ground and first ex- 
cited states are shown in Fig. [j] 

DWBA calculations were carried out with the finite 
range code PTOLEMY [jljj , using the full transition op- 
erator. Seven different optical potentials were studied 
for the 3 He — w O entrance channel. Six came from an 
extensive study of 3 He elastic scattering on s — d shell 
nuclei at 25 MeV Q, with small (<0.5%) adjustments 
in the depths to account for the energy-dependence of 
the real and imaginary volume integrals fi~9[| . One came 



from a global fit |L9|] . The potentials include three differ- 
ent families of discrete ambiguities, characterized by the 
real volume integral, and contain both volume and sur- 
face imaginary forms. In general, the calculations with 
potentials including volume imaginary terms reproduced 
our measured angular distributions slightly better. Even- 
tually, the potentials with the intermediate real volume 
integrals, which were identified as the "physical" family 
in ]18| ], were adopted. The deep potentials predicted a 
forward maximum for the 17 F excited state that varied 
too slowly with angle compared to our measured angular 
distribution. Some of the shallow potentials gave reason- 
able fits to our measured angular distributions at 29.7 
MeV but did a poor job reproducing the 25 MeV data 
Five d — 1( F exit channel potentials were studied. 
Three came from various global fits BOfl, and two came 
from fits to d — 17 O elastic scattering plf. One global 
potential predicted a forward maximum that varied too 
slowly with angle, while the two d — O potentials gave 
very poor fits. The remaining global potentials repro- 
duced the measured angular distributions well and were 
adopted. The single particle orbitals were calculated in 
Woods-Saxon potentials with tq in the range 1.15 — 1.35 
fm and ao in the range 0.55 — 0.75 fm. Over this full 
range, the extracted 17 F ANC's varied by only ±1.5% 
and ±4% for the ground and first excited states, respec- 
tively, demonstrating the insensitivity of the extracted 
ANC's to assumptions about the 17 F wave functions in 
the nuclear interior. In contrast, the more traditional 
spectroscopic factors varied by ±45% and ±19%. 

Normalizing the DWBA calculations to the data and 
accounting for the ANC's for the single particle orbitals 
and the known ANC for 3 He — > d + p [22 2^] provides 
C 2 for 17 F — > w O + p. Fits over several angular ranges, 
from e cm = 2 - 6° to 6 cm = 2-30°, gave ANC's con- 
sistent to within 2%. The final ANC's were determined 
from fits to the forward angle peaks (6 cm = 2 — 9°) to 
minimize the sensitivity to the choice of optical model 
parameters. Table I shows the adopted optical model pa- 
rameter combinations that gave the smallest and largest 
ANC's. It is worth noting that most optical potentials 
that gave poor fits nonetheless gave ANC's that also fell 
within this range. The corresponding fits to the ground 
and first excited state angular distributions are shown in 
Fig. [l]. The fits to the 17 F first excited state near the 
minimum and the weak population that we observe for 
the 17 F | second excited state and | third excited 
state set upper limits on the contributions due to com- 
pound nuclear effects and multi-step reactions at <1%. 
Our final adopted ANC is Cj = 1.08 ± 0.10 fm- 1 for 
the ground state. The uncertainty includes ±4.8% from 
the absolute normalization and angle accuracies, plus the 
statistics of the fits, and ±7.6% associated with the choice 
of optical model parameters and single particle orbital, 
as well as ambiguities in the reaction mechanism. Our 
final adopted ANC is C 2 i/2 = 6490 ± 680 fm" 1 for the 
first excited state. The corresponding contributions to 
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TABLE I. Adopted optical potentials. Sets I and II gave the smallest and largest ANC's for the two transitions, with 
other optical potential combinations giving ANC's in between. The d potentials are specified for the 17 F first excited state. 
Energy-dependent terms were slightly different for the ground state. All energies are in MeV, distances are in fm, and ANC's 
are in fm . 
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its uncertainty are ±5.4% and ±9.0%. 



III. S-FACTORS FOR 16 0(p, 7 ) 17 F 

The relation of the ANC's to the direct capture rate at 
low energies is straightforward The cross section for 
the direct capture reaction A + p — » B + 7 can be written 
as 



A|</* p (r)|0(r)|^ +) (r)>| 



(3) 



where A contains kinematic factors, 1^ is the overlap 
function for B — > A + p, O is the electromagnetic transi- 
tion operator, and is the incident scattering wave. If 
the dominant contribution to the matrix element comes 
from outside the nuclear radius, the overlap function may 
be replaced by 



C 



W_ 77>i+ i/ 2 (2Kr) 



(4) 



where C defines the amplitude of the tail of the radial 
overlap function 1% W is the Whittaker function, r\ is 
the Coulomb parameter for the bound state B = A + p, 
and k is the bound state wave number. For le O(p, j) 17 F, 
the necessary C"s are just the ANC's determined from 
the 16 0( 3 He,d) 17 F transfer reaction studies in Sect. II. 
Thus, the direct capture cross sections are directly pro- 
portional to the squares of these ANC's. In fact, the 
16 0(p, 7) 17 F reaction populating the very weakly bound 
17 F first excited state provides an extreme test of the con- 
nection between the ANC measured in a transfer reaction 
and the S-factor measured in direct capture. The ap- 
proximation of Eq. ([!]) is excellent at large radii, but the 
proximity of the node in the 2si/ 2 wave function makes it 
rather poor near the nuclear surface. In contrast, Eq. (^) 
provides a good description of the 17 F ground state 1^5/2 
wave function even in the vicinity of the nuclear surface. 

Following the prescription outlined above, the S- 
factors for 16 0(p,j) 17 F were calculated with no free pa- 
rameters. The results are shown in Fig. ||. Both El and 
E2 contributions have been included in the calculations, 
but the El components dominate. The capture of pro- 
tons by le O at low energies occurs at very large distances 
r due to the extremely small proton separation energy of 
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FIG. 2. A comparison of the experimental S-factors to 
those determined from the ANC's found in le O( 3 He,d) 17 F. 
The solid data points are from [9], and the open boxes are 
from [10]. The solid lines indicate our calculated S-factors, 
and the dashed lines indicate the ±1<t error bands. Note that 
the experimental ground state S-factor may be contaminated 
by background at energies below 500 keV [25]. 
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17 F |l. Thus, we find that the calculated capture cross 
sections are sensitive neither to the behavior of the over- 
lap functions at small r, nor to the nuclear interaction 
between le O and p in the initial state fllQfl . We find that 
S(0) = 0.40 ± 0.04 keVb for populating the 17 F ground 
state and S(0) = 9.8 ± 1.0 keV-b for populating the first 
excited state. The uncertainties in these calculated zero- 
energy S-factors come almost entirely from those in the 
corresponding ANC's determined above. There is no un- 
certainty associated with ambiguities in an extrapolation 
from higher incident energies to zero energy, and there is 
very little theoretical uncertainty, since the capture re- 
action is almost purely peripheral at very low incident 
energies. In the astrophysical domain, the energy de- 
pendence of the capture cross sections is determined en- 
tirely by the initial Coulomb scattering wave functions 
and the kinematic factors, while their magnitudes are 
fixed by the ANC's. The theoretical uncertainty in the 
S-factors is less than 2% at an energy of 1 MeV. This 
was estimated by repeating the calculation while com- 
pletely neglecting the nuclear interaction in the initial 
state. Hence, the uncertainty in S at small energies is due 
just to the uncertainties in the ANC's measured above. 
However, as the energy increases above 1 MeV, the cal- 
culated S-factors become more sensitive to the behavior 
of the overlap functions at smaller r and to the details 
of the nuclear interaction in the initial state. In that 
case, the simple direct radiative capture model used here 
breaks down, and a microscopic approach including an- 
tisymmetrization is needed. This effect has been studied 
for 7 Be(p : j) 8 B in @. 

Two previous measurements of le O(p,j) 17 F have de- 
termined the capture cross sections to the ground and 
first excited states separately The experimen- 

tal results for the S-factors populating the 17 F ground 
and excited states are also shown in Fig. It is clear 
from Fig. || that the agreement between the experimen- 
tal results and the predictions based on our measured 
ANC's is indeed very good for proton energies below 1 
MeV. At these energies, the 16 0(p, 7 ) 17 F S-factors de- 
rived from the analysis of our le O( 3 He,d) 17 F measure- 
ments agree with the corresponding direct experimental 
results to better than 9%. 

Our calculated S-factors for le O(p,^f) 17 F in Fig. ||are 
very similar to the S-factors calculated for the same re- 
action in H . The energy dependences are virtually iden- 
tical. For both states, we calculate the S-factor to be 
slightly larger than those in which provides us with 
a somewhat better representation of the ground state S- 
factor and a slightly poorer representation for the first 
excited state. It is important to recognize that the pro- 
cedures used in the two calculations are very different, 
even though their final results are quite similar. In 
the 17 F ground and first excited states were assumed to 
be good single particle states outside a closed le O core. 
Electron scattering data were used to specify the density 
distribution of le O, which provided the input for a fold- 
ing model calculation of the low energy p — 16 O potential 



with DDM3Y. The central and spin-orbit terms in the po- 
tential were renormalized separately, for both even and 
odd partial waves, by fitting the 17 F bound state energies 
and comparing to detailed data on low energy p + 16 
elastic scattering. Finally, the direct capture rates were 
calculated with no additional free parameters. This level 
of detail was necessary to reproduce the 16 0(p,j) 17 F S- 
factors at proton energies higher than we consider here. 

The ANC technique is quite different, much simpler, 
and based on much less experimental input. Our mea- 
sured 16 0( 3 He,d) 17 F angular distributions determined 
the ANC's for 17 F -> 16 + p experimentally. These 
specify the amplitudes of the tails of the 17 F — » le O + p 
overlap functions. We then normalized single particle or- 
bitals to the measured ANC's, and used them to calculate 
the corresponding direct radiative capture S-factors. So 
long as one restricts oneself to the low energies typically 
of greatest importance to nuclear astrophysics, the only 
input required by this technique is the experimentally 
measured value of the ANC. In practice, the close agree- 
ment between our calculated S-factors and those in H in- 
dicate that the body of experimental data used to specify 
the p + 16 O potential in |)| ultimately was sufficient to 
determine the 17 F ANC's indirectly. However, the ANC 
approach may also be used to determine S-factors for di- 
rect radiative capture from peripheral transfer reaction 
data in cases, such as radioactive targets, for which much 
less experimental data are available than for le O. 

IV. CONCLUSION 

In conclusion, the 16 0(p,j) 17 F S-factors derived from 
the analysis of our 16 0( 3 He,d) 17 F measurements agree 
with the corresponding direct experimental results to 
better than 9%. This demonstrates the practicality of 
determining accurate S-factors for very low energy di- 
rect capture reactions from measurements of the cor- 
responding asymptotic normalization coefficients in pe- 
ripheral proton transfer reactions. This technique can 
be extended to other systems, including measurements 
with radioactive beams. The production of 8 B in the 
sun via the 7 Be(p, ^f) 8 B reaction is an ideal example. 
While this reaction is relatively unimportant in the pro- 
duction of energy, it provides the only source of high 
energy neutrinos. Hence its rate is crucial to interpret- 
ing measurements from solar neutrino detectors ||. At 
stellar energies this reaction is completely dominated by 
direct capture which occurs at large radii. Indeed, even 
before this demonstration of the accuracy of this indi- 
rect technique, there has been an attempt || to deter- 
mine the 7 Be(p, j) 8 B S-factor from a measurement of 
the S B -> 7 Be + p ANC in the reaction 2 H( 7 Be, 8 B)n. 
But interpretation of that result suffered from significant 
uncertainties in the choice of optical potentials [Q, at 
least in part due to the very low energies involved. The 
8 B — > 7 Be + p ANC can also be measured in ( 7 Be, 8 B) 
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transfer reactions at higher energies with heavier targets, 
where the uncertainties due to the choice of optical po- 
tentials are much reduced. We will report cross sections 
for this reaction using 10 B and 14 N targets in future pub- 
lications. 
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